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ANALYSIS  OF  CONVENTIONAL  AND  REFLECTIVE  BUTLER  MATRICES 
WITH  IMPERFECT  COMPONENTS 


INTRODUCTION 

A  Butler  matrix  that  forms  a  cluster  of  beams  evenly  distributed  in  the  sin-tf  space  is 
not  usually  symmetric  with  respect  to  a  plane  midway  between  the  input  and  output 
ports.  IIo’"  'ver,  by  properly  adjusting  the  phase  shifts  and  interconnections  one  may 
modiP’  conventional  Butler  matrix  to  be  symmetric.  Such  a  matrix  may  also  be  folded 
on  itst-u  on  the  line  of  symmetry,  so  that  the  input  and  output  ports  are  identical.  Such  a 
network  not  only  reduces  the  number  of  components  required;  it  also  becomes  a 
reflection-type  system  in  which  the  feed  positions  are  in  the  plane  of  the  aperture.  The 
synthesis  of  this  network  was  described  previously  [1,2] .  In  this  report,  we  analyze  the 
performance  of  both  conventional  and  reflective  Butler  matrices.  In  particular,  we 
investigate  the  effect  of  reflected  waves  on  the  beam-forming  performance.  In  a  conven¬ 
tional  Butler  matrix,  since  the  input  and  output  ports  are  separate,  the  reflected  waves 
emerging  from  the  input  ports  have  no  effect  on  the  beam -forming  performance.  Multiply 
reflected  waves  may  emerge  from  output  ports;  however,  their  amplitudes  are  generally 
small,  and  their  elfects  are  relatively  insignificant.  In  a  reflective  Butler  matrix,  the 
reflected  waves  accumulate  at  the  innut/output  ports;  hence,  the  aperture  distribution  at 
the  antenna  array  is  significantly  modified,  and  this  may  degrade  the  beam-forming 
performance.  These  effects  are  investigated,  and  computer  simulated  results  are  presented 
together  with  a  listing  of  the  computer  program. 


SCATTERING  MATRIX  OF  A  3-dB  HYBRID  COUPLER 

The  basic  building  block  of  a  Butler  matrix  is  a  3-dB  hybrid  coupler.  For  the  ideal 
hybrid  coupler,  energy  fed  into  any  one  of  the  input  ports  will  be  split  into  two  equal 
components,  one  with  a  phase  shift  of  90^  relative  to  the  other.  However,  practical  hybrid 
couplers  will  in  general  exhibit  amplitude  and  phase  errors  in  their  transfer  coefficients. 
These  amplitude  and  phase  errors  will  affect  the  transfer  coefficients  of  both  reflective  and 
conventional  Butler  matrices  in  the  same  way.  Thai  is,  the  errors  jri  the  overall  network 
input/output  transfer  coefficients  will  be  the  same  for  both  conventional  and  reflective 
networks.  Practical  hybrid  couplers  will  also  have  nonzero  reflection  and  tranrfer  coeffi¬ 
cients  to  the  isolated  port.  For  the  conventional  network,  to  a  first  order,  the  error 
components  due  to  these  effects  will  appear  at  the  network  inputs.  For  the  reflective 
network,  with  its  inputs  and  outputs  sharing  a  single  set  of  ports,  all  error  components 
affect  the  inpul/output  transfer  coefficients. 

Thus,  the  two  types  of  hybrid  coupler  errors  are  forward  and  reverse,  Our  investiga¬ 
tion  v/ill  be  concentrated  on  the  reverse -error  components,  and  we  shall  assume  that,  there 
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is  no  amplitude  or  phase  error  in  tiie  forward-transfer  coefficients  of  the  3-dK  coupler.  The 
-following  analysis  is  based  on  the  assumption  that,  when  an  incident  wave  of  unit 
amplitude  is  applied  to  one  of  the  input  ports,  two  waves  of  amplitude  a  will  emerge 
from  the  two  output  ports,  one  with  a  90°  phase  shift  and  the  other  with  no  phase  shift. 
Similarly,  waves  of  amplitude  P  will  be  reflected  to  the  two  input  ports.  As  shown  in 
Fig.  1(a),  when  an  incident  wave  of  unit  amplitude  is  applied  at  port  12,  reflected  waves 
of  -p  and  -jfi  appear  at  ports  11  and  12  respectively  and  waves  of  -joc  and  oc  appear  at 
ports  21  and  22.  For  conservation  of  energy,  one  has 

2a2  +  2/32  =  1 .  (1) 

The  isolation  factor  is  defined  as  the  power  ratio  of  the  reflected  wave  to  the  incident 
wave.  In  this  case,  the  isolation  is 

/  =  J32.  (2) 


Accordingly,  in  terms  of  the  isolation  factor, 


/0.5  -  1 


If  the  parameters  in  Fig.  'l(b)  are  used,  the  reflected  waves  are  related  to  the  incident 
waves  by  the  matrix  equation 
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where-  a11,a12<a2'i  >  anc*  a22  are  i'icident  waves  and  b11,bl2, b21,  and  b2  2  are  scattered 
waves  at  ports  11,  12,  21,  and  22  respectively. 
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(a)  Unit-amplitude  wave  incident  at  one  of  the 
input  ports  of  a  3-dB  coupler 


(b)  Incident  and  reflected  waves 
on  a  o-aB  coupier 


Fig.  1  —  Transfer  and  reflecti'"  in  four-port  networks 
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Matrix  Eq.  (4)  can  now  be  simplified  to  the  form 
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SCATTERING  AND  TRANSFER  MATRICES  OF  A  BUTLER  NETWORK 

A  Butler  network  can  be  represented  by  a  block  diagram  as  shown  in  Fig.  2,*  Blocks 
in  regions  1  and  3  represent  the  3-dB  couplers  described  in  the  previous  section,  and  a 
phase-shift  transfer  network  is  located  in  region  2.  A  number  of  similar  networks  are 


the  remainder  of  this  report,  a  network  will  be  considered  a  physical  entity  and  a  matrix  a 
mathematical  entity. 
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Equation  (7)  can  now  be  simplified  to 


The  scattering  matrix  in  region  2,  which  is  a  phase-shift  and  transfer  network,  can  be 
represented  as 


Ruj  ^2 
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where  dx ,  d2 ,  cx ,  and  c2  are  vectors  such  that 
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Matrices  Rj  j  and  R22  are  zero,  and  matrices  Rj  2  and  R21  have  identical  elements.  These 
matrices  describe  the  phase  shifts  and  interconnections  from  one  row  of  couplers  to  the 
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next.  Their  elements  depend  on  the  configuration  of  the  Butler  network.  As  an  example, 
the  R  matrix  of  the  4-port  Butler  network  shown  in  Fig.  3  is 


1^1  2  ~  ^21 


10  0  0 
.n 

0  0  e''4  0 

.77 

0  e1*  0  0 


(12) 


[0  0  0  1 


Fig.  3  —  Four-port  Butler  net  work 


Since  we  are  interested  in  the  overall  scattering  matrix  of  this  network,  we  must  first 
convert  the  scattering  matrix  in  each  region  to  i  transfer  matrix,  which  in  turn  can  be 
multiplied  to  form  the  overall  transfer  matrix  of  the  whole  network.  A  transfer  matrix  can 
be  represented  as 


■a 

1 _ 
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1 2 1  j  l22 
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(13) 


where  a5  and  bj  are  the  incident,  and  reflected  waves  at  the  left  hand  ports  and  a2  and  b2 
are  similar  waves  at  the  right  hand  ports. 


It,  can  be  shown  that  a  matrix  T  is  related  to  an  $  matrix  by  the  following  relations 
[3,4] : 


T  -a  a-1  o 

1  1  1  °21  22  3  2  J  1 


T  =  S  S_1 
J  12  °22  -  12’ 


(14a) 
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Since  Sj  2  =  S21 ,  one  may  use  the  simpler  relation  of  Eq.  (16b)  instead  of  Eq.  (16c). 

Elements  of  matrix  S2l  (or  Sj  2)  represent  the  transmitted  waves  at  the  output  ports 
when  a  unit  incident  wave  is  applied  at  any  one  of  the  input  ports.  Therefore,  matrix  S2] 
is  the  transfer  function  of  a  conventional  Butler  network.  Elements  of  matiix  Sj  5  (or  S22) 
represent  the  reflected  waves  at  the  input  ports  when  a  unit  incident  wave  is  applied  at 
any  one  of  the  input  ports.  In  a  reflective  Butler  network  both  the  reflected  waves  and 
transmitted  waves  emerge  from  the  same  set  of  ports.  Therefore,  the  scattering  matrix  of 
such  a  network  is  the  sum  of  matrices  Sj  2  and  Sj  j ,  or 

S  =  Sn+S12.  (17) 


In  deriving  Ibis  relation,  we  have  made  the  assumption  that  the  symmetry  plane  of  a 
reflective  Butler  network  exhibits  an  open  circuit  unity  reflection  coefficient. 


PATTERNS  OF  AN  ARRAY  FED  BY  A  BUTLER  NETWORK 

Figure  4  shows  a  schematic  diagram  of  a  reflective  Butler  network,  which  has  half 
the  components  of  a  conventional  Butler  network.  There  are  n  ports,  since  ports  a j  1 , 
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Jl  a >  ■■  «  ain  416  identical  with  ports  1  >  °2  2>  •••>  a2a 
and  setting  [b2]  =  Ca^  ]  =  0,  this  can  be  represented 


.  Using  previously  developed  notation 

as 
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(18) 
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Jig.  4  Reflective  Butler  network 

J  he  vector  input  of  [a,  ]  can  be  represented,  for  the  ease  of  an  incident  plane  wave 
received  by  a  linear  array,  by 

alk  =  Ak  exP  [JO*  -  1)‘<J  (19, 

where  u  =  2nd  sin  0/A, 
with  A  =  wavelength, 

0  =  angle  of  incidence  from  the  normal  to  the  array,  and 
cl  -  element,  spacing. 
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In  the  subsequent  discussion,  we  shall  assume  that  the  array  has  a  uniform  illumina¬ 
tion  function,  that  is,  that  Ak  =  1.  The  scattering  m8<rix  [Sj  j  ]  +  [S12]  is  computed  as  a 
function  of  isolation  factor  J.  Radiation  patterns  of  the  network-fed  array  are  represented 
by  two  types  of  plot.  One  shows  the  main  beams  formed  by  several  ports  of  the  reflective 
Butler  network,  and  the  other  shows  the  complete  array  pattern  of  one  port  of  the  net¬ 
work,  in  the  range  0  u  <  180°. 


Figure  5  shows  the  array  patterns  of  an  eight-port  reflective  Butler  network.  Figure  5a 
shows  four  of  the  main  beams  for  variation  of  the  isolation  factor  of  the  3-dB  hybrid  from 
10  dB  to  40  dB  Figure  5b  shows  the  array  pattern  when  the  main  beam  is  at  u  =  22. 5C  for 
the  same  range  of  isolation  factor.  Figure  6  shows  the  corresponding  patterns  for  a  16-port 
reflective  netv  ork.  From  these  figures,  it  can  be  seen  that  the  null  filling  level  is  roughly 
equal  to  the  i  .olation  factor  of  the  3-dB  couplers.  That  is,  for  the  case  of  10-dB  isolation, 
the  pattern  i ,  filled  to  a  level  of  about  10  dB  below  its  peak;  and  for  the  case  of  40  dB 
isolation,  the  pattern  is  filled  to  a  level  of  about  40  dB  below  its  peak. 


Tables  1  and  2  show  computed  results  for  eight-port  and  16-port  reflective  Butler 
networks,  respectively.  The  isolation  factors  in  dB  are  listed  in  the  first  column.  The 
transmitted  power  is  the  percentage  of  incident  power,  averaged  over  all  inputs  and  out¬ 
puts,  that  would  emerge  from  the  outputs  for  the  conventional  Butler-network  configura¬ 
tion.  The  remaining  power  emerges  from  the  input  ports.  It  is  seen  that  the  transmitted 
jiower  decreases  as  the  isolation  decreases  and  as  the  number  of  rows  of  couplers  in  the 
network  increases.  For  the  reflective-network  configuration,  the  input  and  output  ports 
arc  combined,  and  the  components  emerging  from  these  ports  are  also  combined.  The  RMS 
amplitude  and  phase  errors  describe  the  effects  of  these  spurious  components  on  the 
combined  outputs  and  are  defined  by 
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where  A b  and  A <p  are  the  RMS  amplitude  and  phase  errors,  respiectively ,  is  an  element 
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<j)k%  is  the  phase  of  and  0^  is  the  phase  of  for  the  ideal  network  with  no  errors. 
The  error  components  increase  with  the  number  of  rows  of  couplers  and  with  decreasing 
isolation. 

A  computer  program  for  carrying  out  these  calculations  is  listed  in  the  appendix. 

In  addition  to  providing  for  imperfect  reverse  parameters  of  the  hybrid  couplers,  the 
program  provides  for  imperfect  forward  parameters  and  for  errors  in  the  interconnecting 
transmission  lines. 


CONCLUSIONS 

An  exact  analysis  procedure  has  been  developed  that  is  applicable  to  both  conven¬ 
tional  and  reflective  Butler  networks  with  imperfect  components.  The  analytical  procedure 
has  been  programmed  for  computation  of  results  for  conventional  and  reflective  Butler 
networks  of  arbitrary  size.  Results  are  presented  for  eight-port  and  16-port  reflective 
networks  using  hybrid  couplers  with  varying  degrees  of  isolation.  The  results  are  given  in 
the  form  of  radiation-pattern  factors  that  would  be  obtained  from  a  linear  antenna  array 
fed  by  the  network  and  also  in  terms  of  the  RMS  phase  and  amplitude  errors  of  the  net¬ 
work  transfer  coefficients, 
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I'  ig.  5a  —  Main-beam  pattern  of  an  eight-port  reflective  butler  network 
isolation  factor  varies  from  10  db  to  40  db 
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I’ig.  6b  —  Array  pattern  of  an  eight-port  reflective  Butler  network 
isolation  factor  varies  from  10  dB  to  <10  clB;  main  beam  at  u  **  22.50 
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Table  1  —  Computed  Statistical  Parameters  for 
Eight-Port  Reflective  Network 


Isolation 

(dB) 

Transmitted 
Power  (percent 
of  incident) 

RMS  Amplitude 
Error 
(percent) 

RMS  Phase 
Error 
(degrees) 

10 

54,85 

30.25 

38.69 

15 

80.78 

19.42 

23.53 

20 

93.20 

12.30 

12.97 

25 

97.77 

7.30 

7.16 

30 

99.29 

4.21 

3.99 

35 

99.77 

2.39 

2.24 

40 

99.93 

1.35 

1.25 

Table  2  —  Computed  Statistical  Parameters  for 
16-Port  Reflective  Network 


Isolation 

(dB) 

Transmitted 
Power  (percent 
of  incident) 

RMS  Amplitude 
Error 
(percent) 

RMS  Phase 
Error 
(degrees) 

10 

44.92 

47.29 

50.58 

15 

76.30 

31.96 

30.95 

20 

91.49 

20.41 

12.66 

25 

97.19 

11.67 

6.91 

30 

99.10 

6.61 

3.84 

35 

99.71 

3.73 

2.16 

40 

99,91 

2.10 

1.22 

Appendix 

COMPUTER  PROGRAM  FOR  ANALYSIS 

This  computer  program  computes  the  coupling  coefficients  from  the  input  ports  to 
the  output  ports  and  the  power  transmitted  and  reflected;  it  also  plots  the  array  radiation 
pattern  if  it  is  desired.  The  type  of  Butler  matrix  analyzed  by  this  program  can  be  either 
a  conventional  or  a  reflective  type  as  described  in  this  report.  For  this  program  three  input 
data  cards  are  required.  The  first  data  card  enters  the  following  fixed-point  (15  format) 
data: 

NPT  —  Number  of  ports  of  the  Butler  matrix  to  be  computed. 

NROW  —  Number  of  rows  of  this  network. 

KLL  —  Absolute  value  of  KLL  represents  the  beam  index  whose  pattern  is  to  be 

plotted.  If  KLL  =  0,  there  is  no  plot.  If  KLL  is  less  than  0,  the  program 
plots  the  array  pattern  and  also  plots  all  main  beams  formed  by  the 
Butler  matrix  network. 

LPR.INT  —  Printout  control.  If  LPRINT  =  0,  the  program  prints  all  detailed  output  at 
each  computation  step. 

The  second  data  card,  which  is  also  in  a  fixed-point  15  format,  specifies  the  number 
of  ports  in  each  basic  coupling  network  in  each  row.  This  implies  that  identical  coupling 
networks  are  used  in  each  row.  However,  coupling  networks  of  different  ports  may  be  used 
in  different  rows. 

The  third  input  data  card,  which  has  a  FI  0.6  floating-point  format,  specifies  the 
coupling  coefficients  of  the  3-dB  coupler  used  as  the  basic  building  block  of  the  Butler 
matrix  network.  These  coefficients  are  read  m  the  sequence  Al,  Bl,  Cl,  Dl.  These 
numbers  are  related  to  the  coupling  coefficient  of  the  3-dB  coupler  by  the  relations  (see 
Fig.  la) 

p.  =  10~<0  05  x  Al), 

P2  =  10~1°-05  X  Bl)( 

ftj  =  io(006  *  Cl), 

and 

«2  =  10~i0.05  x  Dl)_ 
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0001  PROGRAM  RFBNTX 

C  THIS  PROGRAM  FIRST  FIGURES  OUT  BUTLER  MATRIX  CONNECTION  AND  PHASE 

C  ANGLE  .  COMPUTES  THE  TRANSFER  FUNCTION  AND  THEN  PLOT  THE  PATTERN 
C  MATRIX  LIMIT  TO  THE  SIZE  OF  6* 

C  COMPILED  ON  JULY  13.19T6  BY  J.  K.  HSIAO 

C  REVISED  ON  AUGUST  18,1976  BT  J.  X.  HSIAO 

C  ABSOLUTE  VALUE  OF  KLL  REPRESENTS  THE  BEAM  INDEX  WHOSE  PATTERN  IS 

C  TO  BE  PLOTTED 

C  KLL  =  0  NO  PLOT 

C  KLL  CRATER  THAN  0  PLOT  PATTERN  ONLY 

C  KLL  LESS  THAN  0  PLOT  BOTH  PATTERN  AND  MAIN  BEAMS 

C  LLL  =  1 ,  FULL  MATRIX,  LLL  =  0  REFLECTIVE  MATRIX 

C  LPRINT  =0,  PRINT  ALL  DETAILED  OUTPUTS 

C  IF  LPRINT  NOT  EQUAL  0  NO  MATRIX  MULTIPLICATION  RESULT  IS  PRINTED 

C  IF  LPRINT  LT  0  PRINT  ONLY  THE  TRANSFER  FUNCTION 

0002  CONMCN/C*L/PLTAY<SOO) 

0003  C0MM9N/Ct4/A1,A2,B1 ,82 

0004  DIMENSION  NB P ( 1 6 ) , NBK C  1 6  > 

0005  01  PENSION  MC<8.64),PHA<8,64) 

0006  DIMENSION  SI H  32 , 3  2  ) , S L 2 6 3 2 , 3 2  ) , S2 1  (32 , 32  ) , S22 C 32 , 32  ) 

0 00  T  COMPLEX  S11.S12.S21.S22 

0006  CALL  PLOTSCPLT AY, 500,0.) 

0009  NMAX-32 

0010  KC=0 

0011  1  READ  100. NTP.NR9U, KLL, LPRINT, LLL 

0012  IFCNTP.EQ.OGC  TO  2 

0013  3  READ  100,<NBP<I),I=l.NR0H) 

0014  100  FORMAT! 1615) 

0015  READ  101,  A1.A2, 81,32 

0016  101  FORMAT (8F10-6) 

0011  irCKC.CT.O JCALL  OR IG IN <14 . , 0 .  ) 

0018  KC-KC*1 

0019  NRl=NRJU*l 

0020  CALL  NTWK<NTP«NR1,NBP»NBK»HC,PHA) 

0021  ifclu.gt.oigc  to  4 

0022  CALL  HLFNTX'NTP.hKl,  N8 P , NBK , MC , PH A) 

0023  4  CALL  TRFMTXGNNBX.NTP.NRl.NBP, NBK, MC.PMA.Sll, 512,521, 522) 

002  4  LL -  0 

0025  CALL  PRT0UT(NTP,S2I,Sll,LL,NMAX,LTFP, LPRINT) 

0026  tTFP=l 

0021  IF(KLL.EQ.O)GO  10  1 

00  28  Np  A  V  - 1 

0029  CALL  PATERN  (NTP, S2 1 . SU ,KLL ,NPAV, NMAX  ) 

0030  GO  TO  1 

0031  2  CALL  ENOPLT 

0032  END 
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0001 

C 

SUBROUTINE  PK  TOUK  NT  P,  TRFF,TRFB,LL  ,NPX  ,l.TFP,L  PRINT  ) 

LLCT.O  FOR  BLOCK,  ANC  LLIS  1H6  BLOCK  NUMBER 

C 

LL  =0  FOR  OVERAL  TRANSFER  FUNCTION 

0002 

DIMENSION  TRFFCNMX,NMX),TRFaCNMX,NMX) 

0003 

C0NMflN/Ct4/Al,A2, 81,82 

0004 

CnNM0N/C*6/AMPlE32,3?>,AN&L(32,32),ANGTE32)tTRFF2E32,32),  TR(32,32 
C>, AMPAVC  32 ), ANGAVE  32 ), AHXE32  ), ANXE32 ) , AMP RH S C  32 ) , ANGR HSC 32 ) , 

C  SUMR  E 1824) 

0005 

COMPLEX  TRFF,TRF8,TRFF2,  TR.SR 

0006 

KC  =  0 

000  7 

P I  =  3 . 1 .15926536 

0008 

RAC=180,/P1 

0009 

K6  =  6 

0010 

LLL  =0 

0011 

IFEA1.LL.0..0R.LL.GT.0  )LLL---1 

00  12 

1F( LL.LE . 0)G0  TO  1 

0013 

PRINT  101, LL 

0014 

101 

FORMATE//, 20X, 'THIS  IS  THE  TRANSFER  FUNCTION  OF  BLOCK', 15) 

0015 

GO  TO  4 

oou 

1 

PRINT  111 

0017 

PRINT  106 

0018 

106 

FORMATE//, 20X, 'OVERAL  TRANSFER  FUNCTION') 

0019 

2 

IFE A1.GT.0.)G0  TO  3 

002C 

PRINT  119 

0021 

119 

FORMAT(//«10X, 'ZERO  REFLECTION') 

C 

GENERATE  TRANSFER  FUNCTION  FOR  AN  IDEAL  BUTLER  MATRIX 

0022 

3 

PRINT  124.NTP.A1 

0023 

PRINT  117 

0024 

124 

FORMATE/, 20X, 'NUMBER  OF  P OR T S' , 1 5 , 5 X , ' I SOL A T I ON ( DB ) ' , FI 0 . 4 , / / ) 

0025 

CALL  TRFIOLCNTP) 

0026 

IFELTFP.GT.O )G0  TO  4 

002  7 

IFELPRINT  . GT.O)GO  TO  4 

0026 

PRINT  107./(AMPT(I,J), J=l,NTP).I=l,NTP) 

0029 

PRINT  117 

0030 

PR  TNT  10T,(<ANGLE.l,J),J=l,NTP),I  =  l.NTP) 

0031 

PRINT  117 

0032 

4 

IFEAl«LE.0,)K6-2 

003  7 

CO  60  K=1,K6 

0034 

SUM  =  0  . 

0035 

00  15  1=1, NTP 

0036 

15 

SU"RE  I )  =  0  . 

0037 

IFELPRINT  .GT.O)GO  TO  75 

0036 

GO  TO  E7l,T2,73,74,76,7T)K 

00  3  9 

71 

PRINT  102 

0040 

102 

FORMATE//, 20X, "AMPLITUDE  OF  FORWARD  TRANSFER  FUNCTION') 

0C4  1 

PRINT  117 

0042 

117 

FORMATE/) 

0  04  3 

GO  TO  75 

0044 

72 

PRINT  103 

0045 

103 

FORMATE//, 20X, 'PHASE  ANGLE  OF  FORWARD  TRANSFER  FUNCTION') 
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004  ( 

PRINT  117 

0047 

G 3  TO  75 

0046 

73 

PRINT  104 

0  04  9 

104 

F0RNATC//,20X,'A«PLITU0£  OF  REFLECTIVE  TRANSFER  FUNCTION*) 

0040 

PRINT  117 

0041 

GO  10  75 

0042 

74 

PRINT  105 

0043 

105 

FflRMAT(//,20X, 'PHASE  ANGLE  OF  REFLECTIVE  TRANSFER  FUNCTION') 

0044 

PRINT  ! 1 7 

0045 

CO  10  75 

0056 

76 

PRINT  109 

0057 

109 

FORNATC//.20X, 'AMPLITUDE  OF  THE  RESULTANT  TRANSFER  FUNCTION') 

0058 

ORTNT  117 

0054 

GO  TO  75 

0  0  6  C 

77 

PRINT  110 

0061 

110 

FORNATC//.20X, 'PHASE  ANGLE  OF  THE  RESULTANT  TRANSFER  FUNCTION') 

0062 

PRINT  117 

0067 

75 

00  67  1  =  1  , NTP 

0064 

00  70  J=1.NTP 

0065 

GO  TO(61,62,63,64,65,66)K 

0066 

61 

ANGT(J)=CABS<TRFF< I,J)) 

0067 

ANG  T  2  =  ANGT ( J )»*2 

0068 

SUN  =  SUH+ANGT  2 

0064 

SUNR(J)=SUHR(J)4ANGT2 

0070 

GO  TO  TO 

0071 

62 

TFCLPRINT  .GT.O)GO  TO  70 

0072 

ANGT  CJ)=CANG(TRFFCI»J)  )*RAC 

0073 

GO  TO  70 

0 C  74 

63 

ANGT(J)=CABS(TRFfl< I, J) ) 

0075 

ANGT?=ANGT(J )»«2 

0076 

SUN=SUN+ANGT2 

0077 

SUNRCJ)  =  SUNR(J)«-ANGT2 

0078 

GO  TO  70 

00  75 

64 

IF  ( LPRINT  oGT.OGO  TO  70 

0080 

ANGT (J)=C ANG CTRFBCI »J))*RAC 

0081 

GO  TO  70 

0082 

65 

TRCI.J)  =TRFFCI, J)*TRF8Cl,J) 

0033 

ANGT(J)=CABSCTRCI. J)) 

0084 

ANGT  Z  =  ANGT(J )««2 

0085 

SUNR(J)=SUHRC J)+ANGT2 

0036 

SUR  =  SUN*ANGT  2 

0087 

IF ( LLL  «  GT ,0 ) GO  TO  70 

009e 

AMPTCI, J)=CANGTCJ)-ANPT(I,J))/AHPT(1,J) 

0035 

GO  TO  70 

0090 

66 

ANGTC J)=CANGC TRCI.J ))»RAC 

0091 

IFCLLL.GT.OGO  TO  70 

0092 

AG  =  ANGT (J)-ANGLCI ,  J) 

0093 

ANGLE J. 1 )=BG 

0094 

IFC ABS(AG).LE.lflO.)GO  TO  70 

0095 

NS  I GN  =  1 

20 
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o<m 

IF< AG.GT.0.)NSIGN=-1 

0097 

ANGL<I.J)*NSIGN«<  3  60.-ABS<AG>) 

0096 

70 

CONTINUE 

0099 

IFCLPRINT  .GT.03GO  TO  67 

0100 

PRINT  107,<ANGT<J),J=1.NTP) 

0101 

107 

F<3RNAT<10X<  8F10.4) 

010  2 

47 

CONTINUE 

0103 

KN00  =  N<30<K,2  3 

0109 

IFCKNCO.LE.O  )G0  TO  60 

0105 

PRINT  122, SUN 

0106 

122 

F0RNITC//.10X, 'TOTAL  POWER  OUTPUT FIO .5 ) 

OlOT 

PRINT  123<<SUNR<I3< I=1,NTP) 

0106 

123 

F0RN*T<//,1CX, 'POWER  FRON  EACH  P  ORT  '  , /,  <  1  OX,  1  OF  10 . 5)  ) 

0109 

60 

CONTINUE 

0110 

IFCLU.GT.03GO  TO  7 

0111 

IFCLPRINT  . GT • 0 )G0  TO  a 

0112 

DO  50  L  =  l,2 

0113 

GO  TO  <  5 1 « 5?  )L 

0119 

51 

PRINT  120 

0115 

120 

FOR NAT<//#20X,  'ERROR  FUNCT I ON  ' ,  / / ,  20X.  'ANPLI TUD E '  , / ) 

0116 

GO  TO  53 

0117 

52 

PRINT  121 

0118 

121 

F0RNATC//.20X, 'PHASE  ANGLE'./) 

0119 

53 

00  50  1*1. NTP 

0120 

GO  TO  <54.55  >L 

0121 

59 

PRINT  107.  <AHPT< J , I ) • J* 1 .NTP) 

0122 

GO  TO  50 

0123 

55 

PRINT  107.  CANGLCJ.T.  )*J  =  l»NTP) 

0125 

50 

CONTINUE 

0125 

TF<LPRINT  .LT.O)RETUKN 

0  1  26 

PRINT  111 

0  12  l 

ub  ii  1=1./ 

0128 

DO  57  1*1, NTP 

0129 

XFCL.GT.l.AND.I.GT.DGO  10  58 

0130 

•NGS=0. 

0131 

ANP  S  =0 • 

0132 

ANGX-O. 

0133 

ANPX=0. 

0135 

58 

DO  56  J=1,NTP 

0135 

anps^anp$.anpt<j,i  ) 

0136 

ANGS  =  ANGS«-flNGLCJ, I ) 

0137 

IF< ANPTC J,I).GT.ANPX)ANPX*ANPT<J,I) 

0138 

56 

IF<A8SCANGLCJ,I)).G7.ABSCANGX))AN&X=ANGLCJ,1) 

0139 

IFCL.GT.DGO  TO  57 

0150 

*HPAVCI>=ANPS/NTP 

0151 

ANGAVC I) =ANG5/NTP 

0152 

INV  /  T  1  -  A  U  0  V 

0153 

ANX  < I )  =  ANG  X 

0155 

57 

CONI INUfc 

0155 

IF<L.GT.1)G<3  10  59 
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01U 

MINT  117 

014.7 

PRINT  107,  <ANPAV(K),  14=1,  NTP) 

0146 

PRINT  117 

OH1 

PRINT  107,<AN*jAV(K),K  =  l  ,  NT  P  ) 

OISO 

PRINT  117 

0151 

PRINT  10  7,  (ANX<K), K  =  1,NTP) 

015  2 

PRINT  117 

015? 

PRINT  10T,<ANX<K),K=l,NTP) 

0154 

PRINT  117 

0155 

59 

CONTINUE 

0156 

ANRS-ANPS/NTP**2 

0157 

ANGS=ANGS/NTP*«2 

0156 

ANGSST*Q. 

015  6 

ANP  SST  =  0  « 

0160 

03  80  I=-i,NTP 

0161 

ANGSS*0. 

0162 

ANPSS=0- 

0163 

03  81  J  =  1 «  NT  P 

0164 

ANP$S=ANPSS»<  ANPTC  J,  I ) -A  NP  A  VC  I  >  )  **  2 

0165 

ANGSS  ANGSS*(ANGLC J.I J-ANGAVCI ))**2 

0166 

AHPSST=ANPSST+<ANPT(J« I)-ANPS)**2 

0167 

81 

ANGSST=ANGSSH-<ANGL( J,I)-ANGS)**2 

0166 

AKPRNS<n  =  SORT<ANPS$  /NTP) 

0166 

ANGRNS(I)=SQRT<ANGSS  /NTP) 

0170 

80 

CONTINUE 

0171 

PR  t  NT  107,CANPRNSCK).K-1 ,NTP) 

0172 

PRINT  117 

0173 

PRINT  107, ( ANGRNSCK  )  ,K  =  1 »NTP  ) 

0174 

Af1°SS  =  SQRTCAPRSST/NTP»*2) 

0175 

ANGSS=SQRT<ANGSST/NTP9*2) 

0176 

PRINT  117 

0177 

PRINT  107,  ARPS,  ANGS,  ANP X , ANGX , AMP S S ,  ANGSS 

0176 

7 

IE(LC-GT.O)RFTURN 

0175 

IFILPRINT  , NE  «0  )RE  TURN 

oiac 

PRINT  111 

0181 

111 

FOP,  NAT (  1  HI  ) 

0182 

L3=K6/2 

0183 

6 

00  10  t'l,L3 

0184 

GO  TO  <11,12,13)1 

0185 

11 

PRINT  112 

0186 

112 

FOR  N  A  V<  /  /  ,  2  0  X  ,  '  10  E  AL  CASE',//) 

0 1  R  7 

GO  TO  14 

0188 

1? 

PRINT  113 

0188 

113 

FORNAT(//,20X, 'ACTUAL  CASE',//) 

0180 

GO  TO  14 

0191 

13 

PRINT  114 

0192 

114 

FOR  NAT <//,20X,  'DIFFERENCE',//) 

0193 

14 

00  30  1=1, NTP 

0194 

03  30  J  = 1 , NT P 

0195 

SR  =  CNPLXCO.,Q.  ) 
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0106 
0  197 
0198 
0  199 
0200 
0201 
0  20  2 
0203 
0204 
0200 
0206 
0207 
0208 
0209 
0210 
0211 
0212 
0213 
0214 
0215 
0216 
0217 
0218 
0219 
0220 
022  1 
0222 
0223 
0224 
0225 
0226 
02  2  < 
0228 
0229 


00  40  K=l,NTP 
GO  10  C41,42,43)L 

41  SR  =  $«4TftFF(I,K>4C0NJG';TRFF2<J,K)> 

GO  TO  40 

42  SR=SR*TRC I,K )* CON J G( TRFF2 C J,K)) 

GO  TO  40 

4  3  SR  -  SR  t  TRF8 (  I  ,K  )*C0 N J G < T RFF 2 C J * K ) ) 

40  CONTINUE 

ANGL(I,J)=CANG(SR)*RAC 
30  ANPT(I,J)=CA9S(SR) 

00  20  K=1 i 2 
GO  TO  (2 1 , 22 )K 

21  PRINT  115 

115  F0RNATC20X,  'AMPLITUDE',/) 

CO  TO  23 

22  PRINT  116 

116  FCRMATC//.20X, 'PHASE  ANGLE'*/) 

23  00  20  1=1, NTP 
GO  TO  ( 24 , 2 5  )K 

24  PRIN1  107,  <ANPra,J),G=-l,NTP) 

GO  TO  20 

26  PRINT  107,  <ANCUI,J),J  =  1,NTP) 

20  CONTINUE 
10  CONTINUE 

IF<  KC.GT.O  JRETURN 
TFCA1.LE.0.)RETURN 
PRINT  118 

118  FORNMCIHI.IOX, 'REFLECTION  MATRIX  IS  USED',//) 
00  5  1=1, NTP 
00  5  J  =  1 , NT P 

5  TRFF2(I»J)=TRFF(I,J) 

GO  TO  6 
ENO 
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0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 


0010 

001  1 
0012 
0013 
0014 
0015 
0016 
001  7 
0016 
0019 
0020 
0021 
0022 
002  3 
0024 
0025 
0026 
002  7 
002  0 
0029 
0030 
003  1 
003  2 
0033 
0034 
0035 
0036 
0037 
Ou  3  6 
0039 
0040 
004  l 
0042 
0043 
0044 
0045 
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102 

101 


SUBROUTINE  TRFMTX<NM,NN,NRt,NBP,hBK,MC,PHA,Sll,S12,S21,S22) 
DIMENSION  N8P<16),NBKU6) 

DIMENSION  MCCHR1 ,NN) .PHAENRl ,NN) 

CINE  NS  ION  S11<NP,NM),S12(NM,NP),S21CNN,NP),S22(NM,NM) 

COHMON/CA5/YllC32.3?>.T12<32,32).T2102o32).T2Z<32.32> 
COPMON/C46/R11C32,  32  >.R  1202, 32),  R21 132,32) ,R22( 32,32) 

Cl  PENSION  SSlU6,B),SS12(a,6),SS21(8,a),SV>2<8,8) 

OIMENSION  PC  T<  32  ) 

COPPLCX  S l 1 , S' 2,$21fS22,Tll,T12,T21,T22,Rll,R12 ,R21,R22,AR,SS11, 
CSS12,SS21,SS22 
C  FIRST  INDEX  ROW 

C  SECOND  INDEX  COLUPN 

00  10  1=1, NR1 

C  TRANSFER  MATRIX  IN  CONNECTION  REGION 

00  11  L  =  1 , NN 
LL=PCCI,L) 

MCT(LL)=L 

PRINT  102,(MC<I,L),L=1,NN) 

PRINT  102.CPCTCL),  U=1.NN) 

PRINT  101.CPHAU.L  ),L  =  1,NN) 

FORMA. C//.C10X, 815)) 

FORMATE//, (10X,8F10«4)) 

00  20  J= 1 , NN 
00  20  M1,NN 
Tll(,l,iO=CMPLXCO.  ,  0,  ) 

T12CJ,X)=CPPLX<0.,0.) 

T21 ( J.K  )=CPPLX(0. ,0. V 
T22(J,K)=CPPLX(0.,0.) 

IF(PCT(J).NE.K)GO  TO  20 
Tll(J*K)=AR(PHA(l,J)) 

T22(J,K)=C0NJG(T11CJ,K)> 

20  CONTINUE 

PRINT  100,<(Tll(P,N),N=l ,NNj ,  P  =1 ,NN) 

PRINT  100,((T22(M,N),N=i,NN),M=l,NN) 

IFCI.CT.DGO  TO  21 
00  22  J-l.NN 
00  22  K  =  1 , NN 
R11(J,K)=T11CJ,K) 

R12(J,K)=CHPLX<0.,0.) 

R21(J,K)=CMPLX(0.,0.) 

P22(J,K)=T22CJ,K) 

GO  TO  23 

CALL  NT  XML  T( NN  »NN ,  T 1 1 , T 1 2 , T2 1 , T2 2 , R 11 , H 12 , R 2 1 , R2 2 ) 

PRINT  100,(<R1 l(M,N),N=l,NN),H=l,NN) 

PRINT  100,C(R12(P,N),N=1,NN),M=1,NN: 

PRINT  100, C(R21(M, N ) , N  =  1 , NN ) , M= l , NN ) 

PRINT  100,E(R22(N,N),N-1,NH)»N=1,NN) 

.00  ry"nHT\»/,vlCX,orlui4.) 

IFU.EQ.NRDGC  TO  10 

TRANFER  MATRIX  IN  BLOCK  REGION 
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24 
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00  4  £ 
0047 
004  8 
0  0  4  S 

005  0 
0051 
00  5  2 
0053 
0054 
0055 
00  5  £ 
00  5  7 
005  8 
0055 
0060 
0061 
0062 
0063 
0064 
0065 
0066 
0067 
0C  6  8 
0065 

0070 

0071 

0C  7  2 
00  73 
0074 
00  75 
00  7  4 
0077 
0078 
0079 
0080 
0081 
008? 
0083 
0084 
0  C  8  5 
0086 
008  7 
0  0  8  P 
0  C  8  8 
00  9  C 
00  91 
0092 
0093 
0094 
0095 
00  9  6 
OC  9  7 
009  P 
00  9  9 
0100 
0101 
0102 
0103 
0104 


>P=NBP< I )  • _ , 

IFO.LF.DGO  TO  26 
If(NP.tg.N6P(l-l))f.C  TO  27 

26  CALL  9LKC 8,NP,5Sll ,5512 ,5S2 1  ,SS22> 

C  RESET  S  MATRIX 

00  24  J=1,NN 
00  24  K=1,NN 

24  S11CJ.K)  =  CMPLX(0.»0  ,  ) 

S 12C  J,K)=CNPLX(0. ,0.7 
S21CJ,K)=CMPLX(0..C.) 

522(J,K3=CMPLX(0.,0.) 

00  25  J  =  1 (NNtNP 

03  25  JJ=1,NP 
Jl  =  JJ-1 
00  25  KK-1,NP 
XJL  =  KK-  1 

SilCJ*J1 »  J  +K  1 ) - S  S 1 l(JJiKK) 

5126  J*  J1 iJ+Kl >  =  S  S l 2  <  JJ»  KK ) 

S21( J*J1 • J*Kl)sSS2 1( JJ'KK) 

S22CJ*J1 ».UK1)=SS22(JJ,KK) 

25  CONTINUE 

PRINT  100iC(Sll(MiN)tN=l»NN)»M=l*NN3 
PRINT  100,((S12CMfN3,N=liNN)»M=lthN) 

PRINT  100,CCS21(M,N),N=1,NN),M-1,NN) 

PRINT  100iC<S22CN,N)tN=I,NN)*M=l,NN) 

C  INVIRSf  S-MATRIX 

CALL  INVS1(NM,NN,S 12 ) 

CALL  5TTRF(NM,NN,S11,S12,S21,S22,TII,T12,T21,T22) 
PRINT  100.aSll<HkN),N  =  I,NN),H  =  J  ,NN) 

PRINT  1CG»C<S12<«»N)»N  1  ?  N  N  3  *  n  -  1  *  N  n  1 
PRINT  100,((S?1(M,N),N-1,NN),M=1,NN) 

PRINT  100.((S22(N,N),N=1,NN),H=1,NN) 

27  00  50  J-l.NN 
OO  5  0  K -  1 »  NN 

T1 1 ( J,K)  =  S11 (J  ,K) 

T12( J,K):512CJ.IO 
T21 CJ,X)-$21(J»K) 

50  T22  C  J,K) =5226  J  ,K) 

CALL  PTXMLTCNM  tNN,  t  11  ,T1  2  ,  T2  1  ,^7,1111,(12,  R2I.R22) 
PR  I NI  IOOiCCRIIC  M(N)»N-1*NN)*M  =  1»NN) 

°R  I  N I  100«C(R12(M#N),N-1 . NN)tM=ltNN) 

PRINT  1 0  0 , ( <  R  2 ’ <  M , N  3 . N - 1 , NN ) , P - 1 , NN ) 

PRINT  100.CCR22<W,N),N=1,NN)*M=J,NN> 

10  CONTI  NUI 

CAL!  INV 5 1 ( NM, NN.R2P > 

no  40  J - l > nn 
oo  40  y- 1 .NN 
Sl?( JtK)-R?2(J,K) 

S2 1 ( J,K) - R 22 C  J i K ) 

Jll C JtK)  =  CMPLX(0.,0  .  ) 

S22( J,K)^CMPLX(0.,0.) 

03  40  L  - 1 tNN 

5I1U.KJ-SII(J|K)-P.2?CJ»L)*R21CL,K) 
S??(J,tO=S??CJ»K)tRl?CJ,L)*P.2?(L,K3 
40  CONTINUE 

PRINT  100 ,(( SI 1( M.N) ,N- 1 ,NN J ,M-1 ,NN) 

PRINT  100tCCS12CM,N3,N=l,NN),M=l,NN) 

PRINT  100fC(S21(HFN)fN=lcNN)fM~lfNN) 

PRINT  100.CC522(«»N),N=I,NN),M=1,NN) 

RETURN 

END 
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000  1 

000? 

0003 

0004 

0C05 

OOOt 

0007 

oooe 

ooos 


0010 


i 

f 

f 

t 


i 

i 


0011 
0012 
0013 
0  0  14. 
0015 
0016 
0017 
0018 
001  s 
0020 
002  1 
0022 
0023 
0025 
0025 
0026 
002  7 
0028 
0029 
00  30 
0031 
0032 
00  3  3 


0036 
00  35 
00  36 
0037 
00  3  8 
00  35 
0060 
004  1 
0042 
0043 


-  -  '  t  Ppj  t  r  ’  tr- 


SUR  ROUT  I NE  STRF(NN,NN,NR!,NBP,N8K,NC,PH»,$11,S12,S21,S22> 
DIMENSION  N8P(16),NBK<1&) 

DIMENSION  MC(NR1,NN),PHA<NR1  ,NN> 

DIRE  NS  ION  SUCNN,NM),S12CNN,NM),S2 lCNH.NM) ,S22(NH,NN) 
C0O«CK/Ct5/lll(8,8),Tl?(8.8),T21C8t8),T22(8,8),RH<8,8),R12C8.8), 
C  R21<8,8),R?2<8,8>,SPACE<716e> 

OIMENSION  MCK32) 

CIMENSION  SS11(2,23,SS12(2,2),SS21(2,?).SS22C2,2) 

COMPLEX  SU,S12,S21,S22,U1,T12,T21,T22,R11,R12,R21,R22.AR,SS1I. 
CSS1 2  «  SS2  1 , SS22 

CALL  TW0PT(SSll,SS12.SS21tSS22,2> 

C  1ST  INDEX, COLUMN 

C  2ND  INDEX, ROW 

00  10  1=1, NR1 

C  TRANSfER  MATRIX  IN  CONNECTION  REGION 

00  11  L  =  1 , HN 
LL=NC<I,L) 

11  NCT ( LL  >  =L 

00  20  J=1,NN 
00  20  K= 1 » NN 
TU(  J,K)=CMPLX(0.,0.) 

T120 J,K)=CPPLX(0.,0.) 

T21<J,K)=CMPLX(0.,0.) 

T22C J,K)=CMPLX(0. ,0.  3 
IF(MCTCJ).HE.K)GO  TO  20 
T11C J,K)=ARCPHACI, J>) 

T22CJ,K)  =  CON.JG(ril  CJ.tO) 

20  CONTINUE 

IFCI.O T,l)GO  TO  21 
00  22  J=1,NU 
00  22  K  =  1 ,NN 
R11C J|K)=T11CJ,K) 

R12C J,K)=T1?( J,K) 

R2IC J,K)=T21( J  ,K) 

2?  R22(J,K)  =  T22CJ  ,K) 

GO  TO  23 

21  CALL  MTXNLTCNH.NN,  T 1 1 , T  1 2  ,  T2 l * T2 2 , R 1 1 , R  1 2 , R 2 1 , R 2 2 ) 

IF( I. EQ, NR  1)00  TO  10 

C  TRANFER  MATRIX  IN  BLOCK  REGION 

C  RESET  S  MATRIX 

23  NPNBP(A) 

IFC  I.LE .1  ) GO  TO  26 
IFCNP.EU.NBPCI-1))G0  TO  27 
26  DO  24  J=1,NN 

00  24  K=1,NN 

24  Sll C J,K)=CNPLX(0. ,0 . ) 

SI2(J,X)”CMPLX(0.,0.) 

S21CJ,K)=CMP1X(0..0.) 

S22C J,K)=CMPLX(0. , 0.) 

CO  25  J=1,NN,NP 


26 
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f 
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0044 

00  25  J  J  =  1 ,  N  P 

004*) 

J1=JJ-1 

1 

0046 

CO  25  KK=1,NP 

1 

0047 

KI=KK-1 

3c 

0048 

S110J4J1.J«-K1)  =  SS1  1CJJ.KK) 

i 

0  0  4  0 

S12L J*J1 , J+K 1)=SS12( JJ,KK) 

| 

005C 

521 CJ4J1 ,J4K 1 )=SS210JJ,KK) 

f 

004  1 

S22C J+ J1 . J+K 1)=SS22<JJ,KK) 

0052 

25 

CONTINUE 

s 

C 

INVERSE  S-WATRIX 

S" 

005  3 

CALL  INVS2(NH,NN,S  12) 

V- 

0054 

CALL  STTRF(NM,NN,S11,S12,S21,S22.T11,T12,T21,T22) 

r 

0055 

27 

DO  50  J -  1 • NN 

0056 

CO  50  L=1,NN 

s 

005  7 

Tit C J.K)=S11 CJ,K> 

V 

0056 

T12L J,K)=S12CJ,K) 

ft 

0  C  5  S 

T21 (J,K)=S21 <J,K) 

V- 

006C 

50 

T22CJ.K)-S22(J»K) 

i: 

0061 

CALL  MTXHLT(NM,NN,Y11. T l 2 # T2 1 , T2 2, R 1 1 , R 1 2 , R2 l , R22 ) 

0062 

10 

CONTINUE 

t 

0063 

CALL  1NVS2(NH,NN,R22> 

0064 

00  40  J  =  1 »  NN 

| 

0065 

00  40  K-1,NN 

». 

0066 

S12(J,K)=R22(J,K) 

f 

0C67 

S21L J,K)=R22CJ»K) 

j 

0066 

SI  1 C J»K)-C«PLXC0.,O. ) 

0065 

522  C J,K) =CHPLXC0. ,0. ) 

0  C  7  C 

to  40  L - 1 i NN 

00  7  ] 

SlKJ»K)  =  SlKJ»K)-R22CJfL)»R21(L,K) 

0072 

S22CJ.X)=S22CJ.K)4RL2<J,L)»R22(LtK) 

OOtj 

40 

CONTINUE 

0  G  7  4 

RETURN 

0075 

FNO 

000  1 

C 

r 

SUBROUTINE  STTRF(N",NN,SU,S12,S21,S22,T11,T12,T?1,T22) 
THIS  SUBROUTINE  INVERSES  S  PATRIX  AND  STORFS  IN  T 

0002 

t 

01  IF  NS  ION  Sll(NM,N*).S12(NH,NH),S21<NN,NN),S22<NN.NH) 

000  3 

DIMENSION  Tll(NN,NM)«T12(NN#NH)fT21CNW*NN)»T22(NN»NH) 

t 

0004 

COMPLEX  SI  1 , S12.S21 »  S22, TJ 1, T12, T21.T22 

1 

0  005 

OC  30  J- 1 t  NN 

t 

0006 

no  3  0  K  =  1 t  NN 

000  7 

T22( J,K)  =  S1?<J  ,K) 

0008 

T2I( J.K)=CMPLX(0- • 0,  ) 

OOOO 

no  30  L= 1 , NN 

0010 

30 

T21(JfK)=T21(J,K)-S12CJ,L)*SU(L,K) 

001  1 

00  31  J- 1 » NN 

0012 

00  31  K  =  1  *  NN 

001  3 

Tl?(  J,K)  =  C«PLXCO..O.) 

0014 

T 1  1(  J,K)=S?l(  J,r.) 

0015 

CO  31  L  =  1  »  NN 

OOU 

T12(J,K)=Tl?(J,X)*S22<J,L)*S12fL,K) 

00  1  7 

31 

T11(J.K)=T11(J,K)  +522L J,L)*T21(L.K) 

0018 

00  20  K  - 1  »  NN 

0015 

00  20  J  ],NN 

0020 

SI  1 ( J.tO-Tll C J, K) 

002  1 

S12(J.K)-T12(J«K) 

002  2 

<21(J,K)=T21(J,K) 

0  C  2  3 

20 

S2?(J»N)=T22(JjK) 

0024 

Rf  TURN 

002  5 

CNO 
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0001 


0  0  02 
0003 
0004 
0005 
0006 
0007 
0000 
0000 
0010 
0011 
0012 
0013 
0014 
0015 
0016 
0017 
0  016 
0019 
002C 
002  1 
0022 
002  3 
0C  24 
0  0  2  5 
0026 
0  0  2  7 
002  8 
002  9 
00  30 
00  3  1 
00  32 
00  33 
0034 
00  39 
00  36 
003  7 
00  38 
0  0  3  9 
0040 
004  1 
0042 
0043 
004  4 
00  4  5 
0  0  4  6 
004  7 
C  0  4  8 
0  0  4  S 


c 

c 

c 

c 

c 

c 

c 


100 

101 


10 


20 


30 


40 


50 


SUBROUTINE  NTXMLT(NN,NN»Rll,R12.R21,R22,Tll,T12,T21.T22) 

THIS  SUBROUTINE  MULTIPLE  SU8HATRICES  R*T  THEN  STORE  THE  RESULT 
IN  R 
S  =  R  *T 

SI  1 =R11*T 1 1 +R1 2«T2  1 
S12=R11*T12*R12»T22 
S21=R21*T114R22*t?1 
S22=R21*T  12*R22*T22 
0 1 -E NS  ION  TTl(32t32)iTT2(  32*32) 

M MENS  ION  TlUNM,Nr.3,T12CNM,NP)  , T2 t < NK, NN) , T 2 2 < NM , NM ) 

DIMENSION  R  1  KNM,NH),R  1 2 (NM ,NM )  ,R2 1 < NM ,NN > , R 2 2 < NM  , NN  ) 

COMPLEX  T11,T12»T21,T22,R11,R12,R21,R22,TTI,TT2 
PRINT  101 

PRINT  100f((Rll(M,N3*N-ltNN)*M-l»NN) 

PR  TNT  ;00.C(R12<m.N),N  =  1,NN)-H=1.NN1 
PRINT  100,((R2l(M,N3,H=l,NN),P=l,NN) 

PRINT  100,((R22(M,N),N=1,NN),N=1,NN) 

PRINT  100.(CTU(M»N),N  =  1,NN).H  =  1.NN> 

PR  I  NT  100,(<T12(MtN),N  =  l,NN)*«=l,NN) 

PRINT  100,(<T21(N,N),N=ltNN),N=l,NN) 

PRINT  100,C(T22(M,N),N=1,NN),M=1,NN> 

PRINT  101 

FORMATE //,< 1  OX, 8*10. 4)  ) 

FORMATE//,  ' . . . .  ',//) 

00  10  J-l.NN 
00  10  K  =  1 • N  N 
TTIC J,K)=CNPLX<0.,0.) 

TT  2( J,K)=CMPLXCO. , 0. 3 
00  10  L-l.NN 

TTUJ,X)=TT1(J,I04R11(J,L>*T11CL,K)*R12(J,L)*T21CL,K) 

TT2<J,X)=TT2EJ,K)tRll<J,L)*U2Cl»K)4R12CJ,L)*T22<L*K) 

CONTINUE 
00  20  J-l.NN 
CO  2 C  K  =  1,NN 
R 1 1 ( Jt  K ) =TT  1  ( Jt  K  ) 

R12CJ,K)-TT2(J,K) 

00  30  J  =  1 • NN 
00  30  K=1,NN 
Til (J,K)-CMPLX<0. , 0-  ) 

TT2(J,K) =CMPLX<0.,0.) 

00  30  C  - 1 1 NN 

TT1<J,M)=TT1(J,K)4R21CJ,L)*T11CL,K)4R2?(J,L)*T21(L,K) 

TT2(J*K)  =  TT2(J,K)4R2lCJ,L)»I12(L*K)«-R22CJ,L)4T22(L,K) 

CONTINUE 
00  40  J  =  1 »  NN 
00  40  K  =  1 , NN 
P.21  <J,K)=TTI  ( J  ,  K  ) 

R22C J,K) =TT2(J,K) 

00  50  J  =  1 »NN 
UO  50  K  = 1 , NN 
Tll(J.X) =R 11 ( J ,K) 

T12(J,K)=R12(J,K) 

T21( J,K)=R21 L  J  ,  K  ) 

T22C  J,X) -R22  ( J ,K) 

PL  TURN 
END 


i 

i 

f 


% 

•» 

Jt 


1 

I 

1 


1 

n 

'4 


.1  tolli  <UbWi .  liitiilliL  jU'iiiiUi,  ».i  >jliiLTa.£jjikjld k ,«di iK3.tM<4LiwUb*L£: .  1^1 ; Jl  JUi isii-, 


NRL  REPORT  8392 


E 


V 

t? 


i 

tr 

l ‘ 


! 

A 


0001 

SUBROUTINE  PATERN  (NTP . TRFF.  MIFB , KLL ,NPA V , NNX 3 

c 

ABSOLUTE  VALUE  OF  KLL  REPRESENTS  THE  BEAN  INDEX  WHOSE  PATTERN  IS 

c 

TO  BE  PLOTTEO  - 

c 

KLL  =  0  NO  PLOT 

c 

KLL  GRATER  THAN  0  PLOT  PATTERN  ONLY 

c 

KLL  LESS  THAN  0  PLOT  BOTH  PATTERN  AND  RAIN  BEANS 

0002 

COMNCN/Cil/PLTAT(SOO) 

0003 

01NENSI0N  TRFF(NNX»NNX3»TRFRCNNX»NMX3 

0004 

COHNON/CtS/PE AKC64 .  1 00  3  ,  PH  AXE  64  3  , P AV ( 64  3 , P N A  V ( 64  3 , K I  NO l< 64 3 , 

C  KIN02(64J,PEAKDBC 1 00 ) , S P ACE ( 1 3 7 2 ) 

0005 

C0NHCK/C*6/CONTA<40S6),SINTAC4096> 

0004 

COHPLEX  TRFF«TRFB,S 

0007 

ZEX3=10.*AL9G10EX3 

oooe 

PRINT  104 

0005 

104 

FOR  NA  I  (  1H  1 ) 

0010 

PI--3 .1415926536 

001  1 

ATR-PI/18C. 

0012 

KPL0T=IA8S(KLL) 

0013 

NTP  2=NTP / 2 

c 

PLOT  FRAHE 

0014 

VSL=80. 

0015 

NY=YSl 

0016 

XSL-180. 

0017 

NX  =  XSL 

0018 

HN  =  5. 

001S 

ST  =  2  • 

O02C 

XN= 10. 

002  1 

YH  =  5. 

0022 

YS=2  . 

0023 

YSH  =  Y5*-YM 

0024 

NT  A  =20»NTP 

0  CJ  /  5 

NT  A l =  NTA  t 1 

0026 

TAINOPI/NTA 

OC  2  7 

PNOR  =  NTP 

0028 

CALL  PHASANCTAINC.IKA) 

0025 

KL=1 

0030 

1FCKLL.LT.0)KL=2 

00  3  1 

KTAl N=2*NTA» 1 

0032 

00  1  I  =  1 1 NT  P 

0033 

00  1  J  =  1 1 N  TP 

00  34 

1 

TRFFE I, J3  =  TRFP( I, J  3  +  TRFBE I, J3 

0035 

00  20  IL  =  l.Kl. 

0036 

IFCIL.LT  .2300  TO  25 

00  3  7 

CALL  PLOT  CXN+4*  ?0  . i  “3  3 

003  8 

NT  A 1 N  =NT  A 1 

0  0  3  8 

K  S  I G  N  =  1 

004  ' 

■)  r 
t,  j 

CALL  FRANCE  X“*TR*XSL*V  31.  *  ST  $  Hii  t  NX  ?  NY  ) 

004  1 

00  20  K-1.NTP2 

0042 

KX  =  0 

0043 

KFL  A  0  =  0 

SHELTON  AND  HSIAO 


t 


£ 

f 


0044 

KAI  =  1 

0045 

KPCONT  =  1 

00  46 

KRIN0=1 

0047 

KRA=0 

0048 

L£  OGE  =  0 

0049 

05  30  IJ=1.NTA1N 

0  0  5C 

IF (  IL.GT.13G0  TO  23 

0051 

NSIGN=-1 

0052 

IF (  IJ.GF.  .NTA13NSIGN  =  1 

0053 

I=IJ-NTA1*NSIGN 

0054 

GO  TO  24 

0055 

23 

t  =  I  J 

005  6 

24 

11  =  IA8SC I ) 

0057 

11=11-1 

0056 

PAR  =0 . 

0055 

PA  I =0  • 

0060 

IFC IL.LT.23GO  TO  2 1 

0061 

IFC  I.LT.KIN01<K3.5R.I.GT.KIND2CK>3G0  TO  30 

0062 

IFCI  .EQ.KIN01CK331 »  ' 1 

0063 

21 

OO  40  J= 1 • NT P 

0064 

S  =  TRFFC J  *K ) 

0065 

JI=(J-1)*I141 

0066 

JMOOsMOOC  JI • IK A ) 

0067 

IFC JROO.EQ.O )JK00=IKA 

0066 

PAR=CONTAC JHO0)*RE ALC  S  3  - S  INI  AC  J ROD  3  ♦  A I R  AGC  S  )  *NS IGN+PAR 

0064 

PAI=CeNTaC.IP003*AI««G(53+SINTA< JM0034RE»L(S)*NS2GN+PAI 

0070 

40 

continue 

0071 

PAT=PAR4*2+PAI**2 

0072 

PAT=PIT/PKCR 

00  7  3 

IF C 1 L • EC • 2 3GC  TO  22 

00  i  4 

IFCIJ.LE .1300  TO  31 

0075 

IFCMT-PAT  1  3  32.31  .  3  3 

c 

fxarine  if  a  raxihur  is  passed 

0076 

32 

K«I  =  1 

C  0  7  7 

IFCIJ.EQ.2)LE0GE=1 

0078 

IFCKRA.LE.03G5  T5  31 

0079 

IFCPAT1.LE .PEAKCK.KPCCNT  33  GO  T5  34 

0  0  0  c 

KINDI CK3=KRIND 

0081 

kflag=i 

0082 

KPCONT  =KK  +  1 

0083 

34 

KMKK  ♦  1 

0084 

?EAKCK.KK3=PAT1 

0085 

KNA  =  0 

0086 

GO  TC  31 

r 

EXAMINE  IF  A  RINIHUR  IS  PASSED 

0087 

"33 

KRA  =  1 

0088 

IFC  K  R I »LE . 0  3  GO  TO  31 

0089 

35 

KRTND=I-1 

0090 

IFCKFLAG.GT. 0 3K I NO 2C K 3 = I - 1 

0091 

KFLAG=0 

30 


-,=.r  ^ 


I 

i 

i 

A 

1 


;•»  + 
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009  2 

KM  n  0 

0093 

31 

PAT1=PAT 

C 

PLOT  PATTI  RN  P  OR  A  GIVEN  3fcAM 

00  9  9 

IFCK.NE.KPL0T3&0  TO  30 

0095 

IFCXJ.LT. NTA13G0  TO  30 

0096 

22 

Db=2CPAT) 

009  7 

T  =  (  1  .  +  OB7TSL  )  *  Y  M  ♦  S  T 

0096 

IFCY .GT.YSH)Y=Y$M 

0  C  9  9 

IFCY.LT. SY}Y=SY 

010C 

P  =  I  -  1 

0101 

X~  PY  xm/nt a 

0102 

iFcn.EQ.nco  to  3 

0101 

CALL  PL0TCX.Y.2) 

0  105 

GO  TO  30 

0105 

3 

CALL  PL0T(X,Y,3) 

010  6 

30 

CONTINUE 

0107 

IFCIL.GE.23GO  TO  20 

0108 

TFCKMA.GT.03G0  TO  42 

010  9 

IFCKPCONT  .F0.KK3K IN02CK)=NTA1 

one 

TFCKTNO2CK3.LF.0  3K  TN02CK3  =  NTA1 

C 

DELETE  THE  MAIN  LOSE 

0111 

GO  TO  53 

0112 

52 

IFCLEDGE.LE.03GO  TO  43 

0113 

KK=KK+ 1 

0119 

PE  AKC  K»  KK)=P ATI 

0115 

53 

DO  44  1  =  1 f  KK 

0116 

55 

PE8KOBCI)=10.*ALOG10CPEPKCK,I3) 

0117 

PRINT  102.K 

one 

102 

FORMAT  C/7*  20  X  . 'SEAM  INDEX', 15) 

0119 

PRINT  101.CPEAK09C13.I-1.KK) 

0  1  2  C 

K  K  =  K  K  -  1 

0121 

00  53  L=1,KK 

0  12? 

IFfl  .1  T.KPC.ONT  3 GO  TO  53 

0123 

pEAKCK,L)-PEAKCK,Ly13 

0125 

S3 

CONTINUE 

0125 

PMAXCK)=0 

0126 

PSUM  -0. 

0127 

PR  INI  101 .CPFAKCK, I) ,I=1,KK) 

0176 

101 

FOR MA TC/V, Cl  OX, 10E 12.4)3 

0125 

00  50  L=1,KK 

0130 

IFCPEAK(K,l).GT.PMAXCK)>PHAXCK)=PEAKCK,L  ) 

0131 

PSUM  =  PSUM  +  PFAKCK,L  3 

0132 

50 

CONTINUE 

0133 

PAVCK  )  =  PSUM/KK 

0135 

PMAXCK)=2CPRAXCK3) 

0135 

PAV(K)=2CPAVCK3) 

0136 

PRINT  103. PMAXCK), PAVCK) 

0137 

103 

format (//» iox. 'peak'.fio.4.5x.  'average  r»Fio.4) 

0138 

20 

CONTINUE 

0139 

RETURN 

0150 

END 
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0001 

0002 

0003 

0004 

0000 

OOOt 

000? 


0008 
ooos 
0010 
0011 
0012 
00  13 
0014 
0015 

oou 
001  ? 
0018 
001S 


0020 

0021 

OC22 

0023 

0024 

0025 

0026 


002? 
0028 
0  0  2  9 
0030 
00  31 
0032 
0033 
0034 


0035 
0036 
00  37 
00  18 
0039 
0040 
0041 
004? 


SUBROUTINE  HLFM I X C NT P , NR  1 , NB P, NBK . MC , PHA ) 

DIMENSION  N8PC16),NGK<16> 

DIMENSION  MC<NR1,NTP),PHA(NR  1,NTP> 

COHKON?Ci3/MCT(64 ) 

OP1ENSION  ANGC64)  ,ATEMP(64> 

NN=NRl/2 
LL=CNRl*  11/2-NN 

C  LL  =  1  NURSE  R  OF  ROWS  IS  EVEN 

C  LL  =  0  NUMBER  OF  ROWS  IS  000 

CALL  PHASUMCNR1 ,NTP,N8P,MC,PHA,ANG) 

00  10  1=1. NIP 
1 1  =  I 

C  O  11  J= 1 »  NR  1 
KK-NC(J.II) 

IF( J.EQ.NN)KKP=KK 
1 1  =  KK 

11  CONTINUE 

C  FIND  THE  JOINT  POINT  THEN  STORE  IN  MCT  ARRAY 

00  12  J  =  1 . N N 
KKS=MC( J.KK) 

12  KK=KKS 
MCT(KKP)=KKS 

C  AVERAGE  THE  PHASE  ANGLES  FOR  SYMMETRICAL  MATRIX 

00  13  J=1.NN 
J J  =  NR  l-J  +  l 
I HC  =MC ( J  » I  ) 

AVG  =  CPHACJ,  IMO  +  PHACJJ.  I)>/2. 

13  PhA(J,IMO-AVG 
PHACJJ. I)=AVG 

10  CONTINUE 

C  CORRECT  PHASE  ANGLE  OF  THE  MIDDLE  RON  WHEN  THE  NUMBER  OF  ROWS  IS 

C  EVEN 

IFCLL.LE.OIGO  TO  1 
N1  =  NNM 
00  20  1=1, NTP 
II=MC<N1,I> 

TN=MC(N1,MCT(I)) 

ATEMPCII )“PHA(NI,I I  ) 

IFCPHACNl,  IN).GT.ATEMPCH)  )ATEMPCII)=PHA(N1,IN> 

20  CONTINUE 

C  CORRECT  THE  PHASE  ANGLE  6Y  AOC TNG  THE  SAME  EXTRA  PHASE  TO  EACH 

C  PORT  IN  A  BLOCK 

NMP=N0PCN1  ) 

NMB=NBK(N1 ) 

CO  21  1=1. NM8 
I  MB  C I - l )*NMP 
AA  -0  . 

00  22  J-l.NMP 
KK=  IMB  +  J 

A=ATEMP(KK)-PHACN1,KK) 


0043 
OG44 
0045 
0  0  4  £ 
0047 
004  8 
0  0  4S 
0050 


IF( A.GT.AA)AA=A 

22  CONTINUE 
IFCAA.LE.O.l  GO  TO  21 
DO  23  J “ 1 . NMP 
KK=IMBiJ 

23  PHA( N1 ,KK )=PHACM  ,KK  )»AA 
21  CONTINUE 

RE  TURN 

C  CORRECT  THE  PHASE  angles  eor  the  case  when  The  number  of  rows  is 
C  COD 


I 


I 
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0051 

1 

CALL  PHASUM<NR1,NTP#NBP«NC.PHA.ATEMP) 

0052 

08  30  1=1. NTP 

0053 

AA-ANGCD-ATEAPtl) 

0054 

jj* ac  <i.n 

0055 

30 

PHA(t.JJ)=AA 

005  6 

PMACHR1 . I >=AA 

0057 

RETURN 

0056 

END 

0001 

SUBROUTINE  PHASUAO NR I.NTP.NBP.MC.PHA.AS) 

0002 

DIMENSION  NBPC161 

0003 

01  ME  NS  ION  MCCNRl.NTP),PH»CNRl  ,NT P >. ASC NT P 3 

0004 

C 

DIMENSION  LAPC2.64)  <A(64) 

SET  THE  PHASE  SHIFT  OF  THE  BOTTOM  ROW 

0005 

NflOH=NRl-l 

0006 

NN=N8PE NRCW) 

0007 

00  1  J= 1 «NN 

oooe 

LAP62, J)=J 

0005 

1 

AS(  J)=PHA<  1,1) 

0010 

KK  =  NN 

0011 

00  10  I  =  1 • NR  ON 

0012 

IX=NROW-I 

0013 

H  =  II<1 

0014 

NN=  NBP ( 1 1 > 

0015 

IFCII.IE.0)NN=1 

0016 

DO  12  J  =  1 1 NT P 

0017 

tAPCl,J)=LAPC2,J) 

0018 

12 

•  6 J  )  =  AS(J) 

0019 

KN  -  0 

0020 

00  20  L=1.KK 

0021 

LU=t APC 1 «L3 

0022 

DO  21  N  =  1 « NT P 

0023 

TF(ACCIl»N).NE.LL)GO  TO  21 

0024 

JJ  =  N 

0025 

GO  TO  22 

0024 

21 

CONTINUE 

0027 

22 

NMO=MOO( JJ.NN) 

0028 

IF(NMD.EQ.O)NMO-NN 

0029 

DO  30  K= 1 . NN 

0030 

IN0=JJ-NMDH6 

0031 

IF(NN.EQ.1)IN0=JJ 

0032 

KN  =  AN ♦ 1 

0033 

LAp  C  2  > AN )  =  INO 

0034 

30 

A50INlJl  =  #0LLi=HHr.  Cli.LLi 

0035 

20 

CONTINUE 

003  6 

AA  =  AN 

0037 

10 

CONTINUE 

0030 

Rf.  TURN 

0039 

END 
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0001  SUBROUTINE  NTWKC NTP  ,NR 1 «N8P . NBK , «C» PH* > 

C ***»*THI S  SUBROUTINE  FINDS  THE  CONNECTION  OF  A  BUTLER  MATRIX  OR  FFT 
C  GIVEN  THE  NUMBER  OF  ROMS  ANO  THE  NUMBER  OF  PORTS  IN  EACH  BLOCK  IN 
C  EACH  ROW 

C*«**«COMPILCO  BY  J.  K .  HSIAO 

C  ***  **FIR  ST  VERSION  IS  CONPILEO  ON  MAY  3,1976 
C*****NTP,  number  OF  TOTAL  INPUT  PORTS  OR  SAMPLES 
C*****NROW,  NUMBER  0F  ROWS  REQUIRED  TO  PERFORM  THE  TRANSFORMATION 
Cf****NBP  ,  AN  ARRAY  STORES  THE  NUMBER  OF  PORTS  IN  EACH  BLOCK  AT  EACH 
C  ROW.  EACH  BLOCK  IN  A  ROW  HAS  THE  SAME  NUMBER  OF  PORTS 
C«****NBK,  AN  ARRAY  STORES  THE  NUMBER  OF  BLOCKS  IN  EACH  ROW. 

C*****MC,  A  TWO  DIMENSIONAL  ARRAY  STORES  THE  CONNCTIONS  OF  THE  NETWORK. 

C  FIRST  INDEX  Of  THE  ARRAY  REPRESENTS  THE  NUMBER  CF  CURRENT  ROW.  THE 

C  LOCATION  OF  THE  SECOND  INDEX  REPRESENTS  THE  PHYSICAL  LOCATION  OF 

C  THE  PREVIOUS  ROW  WHILE  THE  CONTENTS  OF  IT  IS  THE  CONNECTION  To  THE 

C  CURRENT  ROW 

0002  CINENSION  MC C N R1 , NTP ), PH  AC  NR 1 , NT P ) 

0003  DIMENSION  NFTS <60 ,NBK C 16) ,NBP < 16 3 

C  COMPUTES  THE  NUMBER  OF  PORTS  IN  EACH  BLOCK 

0006  NROW -NR  1-1 

0005  PI=3. 1615926536 

0006  PI2=PI*2. 

OOOT  NBP (NR  1 )  =  1 

0006  NTP2  *NTP/2 

0009  CO  10  1=1, NR1 

0010  10  NBK(I)=NTP/NBP(I) 

C***  *  NFTS  ARRAT  STORES  THE  LOCATION  OF  THE  SAMPLES  IN  EACH  BE  AN ( OR 
C  FREQUENCY  SAMPLE).  THE  STRUCTURE  IS  CHARACTERIZED  BY  TWO  NUMBERS, 

C  NTS, NUMBER  OF  TIME  SAMPLES(OR  INPUT  PORTS)  AND  NFS,  NUMBER  CF 

C  FREQUENCY  SAMPLESCOR  NUNBER  OF  BEAMS).  FOR  EXAMPLE,  NFTSCC3-1.)* 

C  NTSil)  IS  THE  PHYSICAL  LOCATION  OF  THE  FIRST  1IMF  SANPIF  IN  THr 

C  THIRO  FREQUENCY  GKOUP<  OR  OF  THE  THIRD  BEAM), THIS  l'S  REPRESENTED 

C  BY  LMC 

C 

C  SET  THE  INITIAL  NFTS  ARRAY 

0011  00  11  1=1, NTP 

0012  11  MFT  S ( I )  =  1 

C**6*  NT S 1  IS  THE  PREVIOUS  VALUES  OF  THE  NUMBER  OF  TIME  SAMPLESCOR  INPUT 
C  PORTS) 

C***»  NTS2  IS  THE  CURRENT  VALUE 

C***»  NFS1  IS  THE  PREVIOUS  VALUE  OF  THE  NUMBER  OF  FREQUENCY  SAMPLESCOR 
C  BEAMS) 

C*****NFS2  IS  THE  CURRENT  VALUE 

C 

C 

C  SET  THF  INIAL  VALUES  OF  NTS  AND  NFS 

0013  NTS  1 =NTP 

0016  N  F  S 1 = 1 

0015  00  20  1=1, NR1 

C  MM  THF  NUMBER  OF  BLOCKS  OF  THE  CURRENT  ROW 
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C 

NN,  THE  NUMBER  OF  PORTS  IN  EACH  8L 

OC16 

NM=N8K<I) 

0017 

NN -N8P( I ) 

C 

SET  NTS2  AND  NFS2 

OC  1  6 

NTS2=NTS1/NN 

0G15 

NFS2=NTP/NTS2 

c**»* 

THE  ACTUAL  REQUIREO  PHASE  GRADIENT  1 

c 

THE  FIRST  BEAN  IS 

0020 

PAG=PI/NFS2 

c  *  *  * 

THAVAILABLE  PHASE  GRADIENT  FOR  THE  l 

0021 

PSG  PI/NN 

0022 

KA  =0 

0023 

08  30  J=1,NH 

0024 

NOO J=HOQ<J,NFSl) 

0025 

IFCMOOJ.EQ.O )N0DJ=NFS1 

0024 

JJ  =  (J-n/NFSl  +  l 

OC  2  7 

PAGG=PAG*(H00J*2-1 ) 

0028 

PASGO=PSG-PAGG 

00  2  0 

DO  30  K  =  1 «  NN 

0030 

K1=K-1 

0031 

KK=KK4l 

0032 

LHC  =  <MOO.t-I)*NTSl+CK-t  >*NTS2*JJ 

0033 

NCLOC=NFTS<L«C) 

0034 

MC< I»HCL<JC)=KK 

0015 

tF(KK.LE.NTP23GO  TO  31 

0036 

KKI=NTP-KK+1 

0037 

PHA(I«KK)-PHACI«KKI) 

0038 

GO  TO  30 

0038 

31 

IF(PASGO.GT.O.)CO  TO  32 

0  0  4  C 

PHA(I  ,KK',  -  ABS<PASGD3*CNN-K) 

004  1 

GO  TO  30 

0042 

32 

»HA<I,KtO=PASGO*tU 

0043 

30 

CONTINUE 

C 

RECORDING  THE  RRE  GUENC  V  SAMPLE  OR 

0044 

NT  S 1 =NTS2 

0045 

NFS 1=NFS2 

0046 

KK  =  0 

c 

HNS  IS  THE  NUMBER  OF  BLOCKS  WITHIN 

0047 

MNS=HH/NT51 

0048 

DO  40  JM.NFS1 

0045 

JHOD^MOOCJ.MNS) 

0050 

IF( JMCD.EQ.O )J MODERNS 

0051 

JJ^C  J-D/HNS  +  l 

0052 

00  40  K  =  1  «  NTS1 

0053 

KK  =  KM  1 

0054 

40 

NFTSCKK)=-CK-n4NFS14C  JMCO-1  )*NN  +  JJ 

0  0  5  5 

20 

<*uN  i  iriUE 

0056 

RETURN 

0057 

ENO 

KC  W 


IN  EACH  BLOCK  IS 


SAMPLES 


35 


3 


1 

J 


SHELTON  AND  HSIAO 


i 


OOOI 

SUBROUTINE  f RAHECXM,YM,XSL,TSL,SY,HN,NX,NY) 

0002 

COMHONXUi/PLTAYCSOO) 

TJ003 

Y«SY*Y«»SY 

0004 

NLAB=FN*.035 

0006 

HLAS=HLA8+.035 

000  4 

WLAB=4.*KLAB/7. 

0007 

XSCL=XSL/NX 

oooe 

YSCL  =YSL  6NY 

000*7 

0Y= YM/NY 

0010 

Y  =  SY 

0011 

NNY  =  N  Y  <T 

0012 

CALL  PLOTCO. ,SY.3> 

0013 

CALL  PLOT  C  XM«  SY,2  ) 

0014 

CA’.L  PLOUXM,YHSY,2) 

0015 

CALL  PLCT(0„.Y«ST,2) 

0016 

CALL  PLOTCO. ,SY, 2) 

0017 

00  10  1=1,2 

0018 

Y=SY 

0019 

IFC I. 07.1)00  TO  12 

002C 

X1  =  0. 

0021 

X2  =-  .2 

0022 

X3  =-. 1 

0023 

00  TO  13 

0024 

12 

X1*XM 

0025 

X2=XM*.2 

0026 

X3=XN«.l 

002? 

13 

00  10  J=1,NNY 

0028 

CALL  PLOT C  X 1 »Y»  33 

0029 

MOO Y=P00CJ~ 1,10) 

0030 

IFCN0CY.NE.03GO  To  11 

0031 

CALL  PL0TCX2*Y,2) 

0032 

XFCI.GT.13GO  TO  10 

00  3  3 

A=TSCL4(J-l-NY) 

0034 

CALL  hUMBER(-6.5*MLA8,Y-HLAe/2.,HLAB ,A,0. ,4HF3.0> 

0035 

GO  TO  10 

0036 

11 

CALL  PLOTCX3.Y,2) 

0037 

10 

Y  =  T  +OY 

0038 

0  X=  X  M/NX 

0035 

NXX=HX*1 

0040 

00  20  1=1,2 

0041 

X  =  0  . 

0042 

IK  I.GT.l  )G0  TO  22 

0043 

Y1  =  SY 

0044 

Y2=Yl-.2 

0045 

YJsYl-.l 

004  6 

GO  TO  23 

0  0  4  1 

22 

¥  i  =  Y  n  3  “ 

0048 

Y2=Yl+.2 

0049 

Y3=Yl+.l 

0C5C 

23 

00  20  X=l,NXX 

0051 

KK=K-1 

0052 

CALL  PLOT C  X* Y , , 3' 

0053 

MO0x=M30CXK,10 ) 

0054 

IFCMOOX.Ne.07GO  TO  21 

0055 

CALL  PLOT ( X, Y2  »  2 ) 

0056 

1FCI.GT.DG0  TO  20 

0057 

A  =  K  X  *  X  SCL 

0058 

CALL  NUM8ER(X--2.5*WLAB,  S  T-ML  AB* 3 . 0  ,HL A B ♦  A , 0  .  ,4HF3.0> 

0059 

GO  TO  20 

0060 

21 

CALL  PLOT  C  X, Y3 , 23 

0061 

20 

x=x«-ox 

f 


3 


3 


i 
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ft 

,a_'. 


0062 

CALL  SYM80L6 

CREFS.0.,22) 

0063 

35  CALL  SVNBOU 

*90. , 18) 

0065 

32  CALL  PLOT 6  0, 

0065 

END 

XH-17.5*WLAB,-5.*HLA8»SY,HLA$,22HPARANETER  U  IN  DEG 
*WLA6,YM/2. ♦ST-15. «  W..  AB  , KL  A$  ,  18H  ARR  AY  PATTERN  <06>» 

.3) 


0001 


0002 
0  0  0  2 
OU05 
0005 
0006 
ooo  1 
0008 
0  0  0  5 
OOlC 
0011 
0012 
0013 
00H 
0018 
0016 
001  7 

ooie 

0019 
0  02  0 

002  1 
0022 
0023 
0025 
0025 
002  6 
0027 
0028 

0029 

0030 

0031 

0032 

0033 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


10 

15 

18 

23 

26 

27 

28 
35 


C 


39 


2 


C 

3  1 

32 

33 
C 


SUBROUTINE  S IMCX6 I S, OR IG.NN. N AT , ACT, ANS.LK ) 

IOENT  NUN8ER  -  FIO02R0O 

TITLE  -  COMPLEX  MATRIX  INVERSION.  SOLUTION  OF  LINEAR  EQUATIONS 
IOENT  NAME  -  Fl-NRL-SIMCX 
LANGUAGE  -  FORTRAN 
COMPUTER  -  CDC-3800 

CONTRIBUTOR  -  JANET  P.  MASON »  CODE  7813,  RESEARCH  COMPUTATION 
CENTER,  MIS  0 IV I S  I  ON 

ORGANIZATION  -  NRL  -  NAVAL  RESEARCH  LABORATORT  -  WASHINGTON,  t.C. 

20390 

DATE  -  16  OECEMBER  1970 

PURPOSE  -  TO  SOLVE  THE  CONPLFX  MATRIX  EQUATION  AX=B  WHERE  A  IS  A 
SQUARE  COEFFICIENT  MATRIX  ANO  8  IS  A  MATRIX  OF  CONSTANT 
VECTORS.  THE  DETERMINANT  AND  INVERSE  OF  A  ARE  ALSO 
CBTAINEO. 

COMPLEX  SUM,  MAT.ORTG,  ANS,B0,B2  , 85, 86, 88, 810,  P,  11, 813, 815,  CC,CC  2,81 
EQU»VALENCE682*C),6CC,CX),(CC2,CX2> 

DIMENSION  HATCMCT, 1 > . OR IGC NN » 1 ) , ANS6 MCT > , C62 > , CX6 2  ) , CX2C2) 
FJRMAT61X,  2E12.6) 

FORMT  6  25H  THIS  MATRIX  IS  SINGULAR/) 

FOR  MAT 6  1H0 •  '  VALUE  OF  DETERMINANT  IS  " s 2E 1 2 . 6, / / ) 

F0RMATClX,2r 12.6, 5X.2E 12.6) 

FOR  HAT68X,  'ORIGINAL  C ON  ST AN T S * , 2 1 X , 'DE R I V E 0  CONSTANTS*/) 
F0RMATC1HX.6X, 'THE  INVERSE  (BY  COLUMNS)') 

FORMATUHO) 

E0RMATC1H1.6X, 'VALUES  OF  THE  UNKNOWNS') 

FORHATCSX ,  'IDENTITY  MATRIX') 

E3=6~l. 0,0-0) 

85=60.0,0.0) 

811=61.0,0.0) 

ICT  =MC  T 
JSINCrMCT 
HT  =  ACT ♦ 1 


nct=mct+mct 

PU1  ORIGINAL  MATRIX  INTO  MAT 

IF6IS.EQ.0)G0  TO  39 

Kl'HIHlS 

NC  T  =  ICT 

DO  2  J  =  1 , 1 CT 

DO  2  1=1, MCT 

MAT6I.  J)--0RTG6I,  J) 

CONTINUE 

IF6  IS.NE.O )G0  TO  30 

PUT  IDENTITY  MATRIX  INTO  RIGHT  HALF  OF  HAT 

00  32  J=HT , NCT 

DO  32  1=1, MCT 

KAT6I, J)=B5 

DO  33  J=l,MCT 

MAT6 J, JHMCT)=83 

FORM  TRIONGULARIZEO  MATRIX 


5 


I 


i 


a 

A 


j 

1 
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0034 

0035 

0036 

0037 

0036 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0046 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

Cnift 
ww 

0061 

006? 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

007C 

0071 

0072 

0073 

0074 

0075 

0076 
0077 
0078 
0  0  79 


30  JCT»NCT-1  _ =: _ •  • — . 

00  3  J - 1 »  J CT 
KK=J»1 
0070  25 

24  00  4  K=KK»NCT 

B8-MKK,  JVNATCJ,  J) 

00  5  L  =  J  tHf  T 
810=B8*HAT<J,L> 

5  NATCK»L)«MT<K,L)-eiQ 

4  CONTINUE 

C  VALUE  OF  DETER* IN ANT 

3  811 =Bll*NATCJ»J) 

611»=811*NATCNCT,NCT) 

L0W=-NCT 

N0=-1 

C  TO  00  ONE  OR  FORE  SACK  SOLUTIONS 

00  6  N  INC  «NT  »NCT 
JFIN=NCT 
IX  =  0 

C  SACK  SOLUTION 

DO  6  INN»L0K*NO 
MIA6SCINN) 

BO--MT(N.NINC) 

B2=NAICN,K) 

B4*<0. 0,0.0) 

IFCIX)  7 , 2  2  «  7 
22  IX=IX+1 
GOTO  fl 
7  NO  2  =  - JEI N 

CO  9  INN -L  urn  •  K  02 
N  =  I AB  S  < INN ) 

9  B4=y4*NAT(N«N)*NATCN»NINO 
00  --B0-B4 
!f  IN  =  JFIN-l 

6  IFCCU).EQ.0.0.ANO.C(2).EQ.0.0)G0  TO  13 
29  NAT(P,MlNC)=B0/82 

ANS<N)=60/B2 
6  CONTINUE 

00  40  J=HT,NCT 
JJ= J-NCY 
00  40  1=1, NCT 
40  ORIGCI>J.I>-RAT(I,J) 

1FCLK.GT.03RFTURN 
IFUS.EQ.OGO  TO  34 
GO  TO  41 

C  CHECK  FOR  SINGULARITY  AND  TO  SEE  IF  FIRST  TERN 
2  5  JV--J 

CC=NATCJ»J) 

IF( CX(1 > . NT .0.0. OR.  CX(2 3, NF, 0.0)00  TO  12 
11  IF< JV.NE. JSING1G0  TO  14 
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0  0  8  C 

1  3 

PR  I  NT  1 5 

0081 

PRINT  100,  J,(MAI(K,J),K  =  1,MCT> 

00  8  2 

100 

FOKMAT(10X,I5,8F1G.4) 

0083 

PRINT  21,<<MAT(I,  J),I=1,NCT>,  J=1,MCT) 

0084 

RF TURN 

0085 

14 

jv-jvh 

0086 

CC?  =  MT(JV,J) 

0087 

TF(CX2(l),Fa.O.O.ANO.CX2<2T.EC.O.O)GO  TO  11 

0088 

16 

00  17  J J  r  J, NCT 

0  0  8  5 

86=  MA  T ( J , JJ) 

0090 

NAT  (J,JJ)=MAT(JV,JJ) 

0091 

17 

MATE JV, JJ)=86 

0C92 

811 =-811 

009  3 

12 

CONTINUE 

0095 

GOTO  24 

c 

PRINT  SUBSTITUTIONS  BACK  INTO  ORIGINAL  NATRIX 

0095 

45 

00  20  NN V  =  1 *  IS 

0090 

PRINT  27 

0097 

44 

PR  I  NT  2  3 

0098 

00  20  LL  =  1  * N CT 

0099 

B13=(0. 0,0.0) 

0100 

00  19  NM  =  1 » HCT 

0101 

19 

B13=ORIGCLL.NM)*MAT(MM»MCT*NNV)+B13 

0102 

B15=-0RIG(LL,NCr+NNV) 

0103 

PRINT21, B15, B13 

0104 

20 

CONTINUE 

0105 

RETURN 

c 

PRINT  TITLE  -  THE  INVERSE 

0100 

34 

PRINT  26 

0107 

GO  TO  43 

c 

PRINT  TITLE  -  VALUES  OF  UNKNOWNS 

0108 

41 

PRINT  28 

C  1  0  7 

43 

CO  38  J  J  -  lili  NCT 

one 

PRINT  27 

0111 

00  38  11=1, «CT 

0112 

38 

PRINT  10,  NATCII, JJ) 

c 

PRINT  VALUE  OF  DETERMINANT 

0113 

PRINT  1C, 811 

0114 

IF( IS.NE .0)00  TO  45 

c 

PRINT  IDENTITY  MATRIX 

0115 

PRINT  35 

0116 

00  36  K  =  1 »  MC  T 

0117 

PRINT  27 

one 

00  36  1=1, MCT 

0115 

SUW=C0.O ,0.0 ) 

012  0 

00  37  J=1,MCT 

0121 

37 

SU»1  =  0R1&(K,  J  >*NAT(  J,HCTH)  +  SUM 

0 12  2 

36 

PRINT  10, SUM 

0123 

RETURN 

0124 

ENO 

1 

3 

j 

j 


i 

i 


f 


1 

i 

•i 


39 
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0001  SUBROUTINE  BLKCNH, NN, Sit t S12. $21 • S22 ) 

0002  DIMENSION  NBPC  16) ,NBK< 16) 

0003  DIMENSION  HC (8 , 1 6 ) ,PHA<8 ,1 6  ) 

000*  DIMENSION  S11(NH,NM),S12CNK,NH).S21(NM,NM),S.?2CNR,  NM) 

0005  CCNH8N/Ct5/T11C8,8).T12(8,  8),T2H8,8),T2  2(8.8),R11(8,8),R12(8,8) 

C  R2U8,8),R22(8,8),SPACE<7168) 

0006  COMPLEX  S11.S12.S21.S22 

0007  COMPLEX  TU.T12.T21.T22.fUl.R12<R2l.R2? 

0006  IF'NN.GT. 2)GO  TO  1 

0005  CALL  IW0PT(Sll,S12,S2l,S22.NM) 

0010  RETURN 

0011  1  11=0 

0012  N2  =  NN 

0013  3  N2=N  2/2 

001*  IF( N2 .LE . 0 )G 0  TO  2 

0015  11=11*1 

0016  GO  TO  3 

0017  2  00  10  1*1.11 

0018  10  NBPCI)=2 

001S  11=11*1 

0020  CALL  NTHKCNN.II.N8P.NBK.MC.PHA) 

0021  CALL  STRF(NM.NN.II.NBP.NBK«HC.PHA.S11.S12»S21.S22) 

0022  RETURN 

0023  END 


0001  SUBROUTINE  I NVS 1 CNM , NN  ,  $12  ) 

0  002  COMM  ON/C 1 5 /A 1( 32), T( 32, 6*). SPACE C *032) 

0002  OIMENSION  S12CNF.NM) 

0006  COMPLEX  A1,T,S12 

0005  CALL  SI«CX(0,S12,NM,T,NN,A1,1> 

0006  RETURN 

0007  END 


0001  SUBROUTINE  INVS 2<NM ,NN , S 12 ) 

0002  CONMON/C*5/AIC6),T(8,16),SPACEC7920) 

0003  OIMENSION  SI2CNN.NN) 

0006  COMPLEX  A  L , T , S 1 2 

0005  CALL  SIMtX10,S12,NM,l,NN»Al,l> 

0006  RETURN 

0007  END 
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0006 

0007 

oooe 

0000 

001C 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

ocie 

0010 

0020 

0021 

0022 

0023 

"0024 

0025 

0026 

0027 

0028 

0026 

003C 

0031 

0032 

0033 

0034 

0035 

0036 

0037 


1 


2 


SUBROUTINE  TUOPT  (  SU  ,  SI  2  ,  S21 ,  S22  ,  H) 

OINENSION  S11<N,A),S12<A,N>,S21<N,A>,$22<M,N> 
COMPLEX  Sll p  S12 > S2 1 >  S22 
COMMON /CI4/A»BtC»D 
8CO=B+C*0 


IFOfiCO.CT.O. )GO  TO  1 

AR=*0.44C-A*.05) 

IE(A.LE.0.)AR=0. 

I1=SQRT(.5-AR*AR) 

A2»SQRT<,5-AR*AR) 

81  =  AR 


82  =AR 


GO  TO  2 

81*10. **(-A*. 05) 
IftA.LE.O. )B1=0. 

82=10. **<-B*. 05) 

IFCB.LE.O. )82*0 • 

Al=10.**(-C4.05) 

A2=10.4*(-D«.05) 

511<l,l)=81»CMPLXC0.,-l.) 

511<2i2)*Bl*CMPLX(0.,-l.) 

S22Cl,l)«B14CHPLX(0..-l.) 

S22C2#2)-81*CMPLXCO.*-l.) 

511C1,2)=82*CMPLX(-1.,0.) 

S1102«1)=B2*CMPLXO-1.,0.) 

S22C1»2)=82*CMPLXC-1.,0. ) 

S22C2, )  )  =  82*CNPLXC-1 .,0. ) 

5 120  1 ,1)=A1»CMPLX( 1. ,0.) 

S12<2,2)=A1*CMPIXC1..0.) 

52t(l,l)=Al*CMPLXCl..C.) 

S21.  C2.2)  =  A1*CMPLX(  l .  ,o.  ) 

S12Ci»2)=A2*CnPLKtu.«-l.) 

S12C2ii)=A2*CMPLXC0.,-l.) 

S21Clt2)=A2*CMPLX(0.»-l.  ) 

S21C2,l)=A2*CMPLX(Q.«-l.) 

RETURN 

END 


! 


1 

1 


4 


l 

m 
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0001 
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0006 

0007 

0008 

0009 

0010 

0011 
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0013 
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0001 

0002 

0003 

0004 

0005 

CGuu 

0007 

0008 

0009 
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0011 
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0002 
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0004 
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0013 


SUBROUTINE  TRT(MTP,TRFF,LL,NNXl 

CONNON/C86/APNT(32»32), ANGL( 32.32 3.ANGT <323, TR 10(32, 32), $P ACE (4064 
*> 

DIME  NS  ION  TRFF(KMX , KMX ) 

COMPLEX  TRFF.TRIQ 
IF(LL.GT.O)GO  TO  I 
09  10  1*1 , NTP 
00  10  J*1 *  NTP 
10  TRIOCI»J)*TRFF<I,J) 

RETURN 

1  09  20  I=1»NTP 

09  20  J  =  1 »NTP 
20  TRFF(I,J)*TRIO(I,J) 

RETURN 

END 


f 

i 

1 

1 

1 

I 


SUBROUTINE  TRFTOL(NTP)  1 

COMNOK/Ct6/AMPT(32,32),  ANGL(32 , 32 > , ANGT( 32 ) , TRFF2 ( 32 , 3 2 > ,  TR (32,32  .■ 

O,SUMR(2016)  s 

COMPLEX  TRFF2.TR  t 

PI  =  3 . 1415926536  ! 

PI2=P1*2.  \ 

"Ti-tan  inr  i 

l\  I  M™  6UW*#  *  * 

A=SQRT(1./NTP3  J 

P=-P I/NT  P  i 

00  10  1*1, NTP 

PPi(I-l)*P  i 

PR=P*( I-.5 )*2  « 

00  10  J= 1 , NT  P  ; 

PP=*MflO(PP ,P 12 )  I 

RE=A»COSCPP) 

RI=A»SIN(PP)  ’ 

TRFF2CI, J)=CMPLX(RE,RI> 

AMPT(I,J)=A 
ANGLd,  Jl*PP*RTA 

10  PP=PP*PR  * 

RETURN  3 

END  * 


SUBROUTINE  PHASAN  (TAINC.I) 
COMMON/C 46/C0NTA( 40 96), SINT A (4096) 
PI-3.1415926536 
PI2  =PI *2 . 

TA-O. 

1  =  0 

1  1=1*1 

CONTA(I)=COS(TG) 

5IWTA(I)=SIN(TA> 

ta=t  a*tainc 

IF(TA.Gfc.PI2)RETURN 

GO  TO  1 

END 
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0001 
0002 
0003 
0004 
000  5 
0006 
000  7 
0008 


CONPLtX  FUNCTION  AR(AUG) 
ANP*1. 

AG- AUG 

RE=ANP*COS<AG> 

RI*A«P*SINCAG) 

AR=C«PLXCRE,RI) 

RETURN 

CNO 


0001 

0002 

0003 

0004 

0005 

0006 

0007 


FUNCTION  CANGCSR) 
CON  PI.  EX  SR 
A>  -RE  4L( SR  ) 

A2 -AIHAGCSR ) 

CANG=ATAN2CA2,A1) 

RETURN 

END 


